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ABSTRACT 

Published galaxy power spectra from the 2dFGRS and SDSS are not in good agree- 
ment. We revisit this issue by analyzing both the 2dFGRS and SDSS DR5 catalogues 
using essentially identical techniques. We confirm that the 2dFGRS exhibits relatively 
more large scale power than the SDSS, or, equivalently, SDSS has more small scale 
power. We demonstrate that this difference is due to the r-band selected SDSS cata- 
logue being dominated by more strongly clustered red galaxies, which have a stronger 
scale dependent bias. The power spectra of galaxies of the same rest frame colours from 
the two surveys match well. If not accounted for, the difference between the SDSS and 
2dFGRS power spectra causes a bias in the obtained constraints on cosmological pa- 
rameters which is larger than the uncert ainty with which t hey are determined. We also 
found that the correction developed by ICole et alj ([2005) to model the distortion in 
the shape of the power spectrum due to non-linear evolution and scale dependent bias 
is not able to reconcile the constraints obtained from the 2dFGRS and SDSS power 
spectra. Intriguingly, the model is able to describe the differences between the 2dFGRS 
and the much more strongly clustered LRG sample, which exhibits greater nonlinear- 
ities. This shows that more work is needed to understand the relation between the 
galaxy power spectrum and the linear perturbation theory prediction for the power 
spectrum of matter fluctuations. It is therefore important to accurately model these 
effects to get precise estimates of cosmological parameters from these power spectra 
and from future galaxy surveys like Pan-STARRS, or the Dark Energy Survey, which 
will use selection criteria similar to the one of SDSS. 
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1 INTRODUCTION 

In the last decade, the advent of new precise cosmological 
observations coming mainly from measurements of fluctu- 
ations in the temperature of the cosmic microwave back- 
ground radiation (CMB) and the large scale structure of 
the Universe (LSS), have shown a dramatic improvement. 
This has revolutionized our knowledge of the values of the 
basic cosmological para meters which are constrained to an 
accur acy of around 10% jSanchez et alj|200d : ISpergel et alj 
120071 ). 

In particular, the measurements of large scale galaxy 
clustering place important constraints on cosmological pa- 
rameters that complement those from the analysis of fluc- 
tuations in the CMB. Measurements of the galaxy power 
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spectrum constrain the parameter combinations Qmh and 
f2b/fim- The constraint on fi m ft is particularly important as, 
for instance, it breaks a degeneracy in the CMB data and 
allow s accurate determination of fi m (e.g. Efst athiou et all 
l2002h . 

In view of these rapid improvements in the amount and 
quality of the observations the control of the systematic ef- 
fects introduced in the analysis pipeline becomes increas- 
ingly important. New experiments and surveys are being 
planned which will push the level of precision of the con- 
straints on cosmological parameters even further. Thus it is 
necessary to establish how robust the constraints really are 
with respect to the details of the hypothesis implemented to 
establish the link between theory and observations. 

The tension found between the constraints coming from 
the published analysi s of the two - degre e field galaxy redshift 
survey (2dFGRS) of ICole et ail l|2005l ) and the Sloan Dig- 
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Figure 1. Equal area, all sky, AitofT projections of all the galaxies accepted by our SDSS (upper) and 2dFGRS (lower) masks used in 
our analysis. 



ital Sky Survey (SDSS) of lTegmark et alj (|2004h might be 
an indication of p ossible systematic effects. In figure I6a of 
ICole et all (|2005l ). which compares the published estimates 
of the galaxy power spectra, there is evidence of more large 
scale power in the 2dFGRS than in SDSS. This folds through 
and results in t he headline valu e of £lmh = 0.168 ± 0.016 
from 2dFGRS jCole et al.1 120051 ) being lower than that of 
iTegmark et al, | l|2004 ) SDSS, fi m /i = 0.213 ± 0.023. Here, 
one should be cautious as di fferent priors have b een as- 
sumed, but in the analysis of ISanchez et al.l l|2006h . which 
treats each dataset on an equal footing and separately com- 
bines each with CMB data, one sees in their figure 18 that 
the SDSS prefers a substantially higher value of f2 m than 
does the 2dFGRS. In fact, while the 2dFGRS estimate is in 
good agreement with that from the CMB data alone (it sig- 
nificantly tightens the constraint without shifting the best 
fitting value) the SDSS data pull fi m to higher values than 
preferred by the CMB. 



The 2dFGRS and SDSS are the largest redshift surveys 
available and give the most detailed description of the large 
scale structure of the Universe as traced by galaxies. Hence 
to be sure that systematic errors are not biasing the parame- 
ters it is very useful to see each data set subjected to analysis 
by a variety of algorithm s and codes, as happened following 
the IPercival et al.1 ll200ll ) 2dFGRS power-spectrum analysis 
jTegmark et al.ll200ll ). 



iPercival et al.l |2007aT ) computed the power spectrum of 
a combined Main galaxy and Luminous Red Galaxy (LRG) 
sample drawn from SDSS Data Relea se 5 using a similar 
technique to that of ICole et al. I i|2005l ) and found an even 
higher value of the matter density i} m h = 0.23 ± 0.01, sug- 
gesting that the inclusion of the LRGs in the analysis in- 
crease the discrepancy between 2dFGRS and SDSS. Their 
results also show that if the analysis is restricted to large 
scales (0.01 < k < 0.06 ft Mpc -1 ) t he data favours a lower 
matter density Q m h = 0.16 ± 0.03. IPercival et"afl (|2007al ) 
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suggested that these differences could be explained by scale- 
dependent galaxy bias on large scales, but found no signifi- 
cant evidence of it. 

Here we seek to investigate whether these differences 
are as a result of: a) larger than expected cosmic variance, 

b) systematics due to differences in the analysis technique 

c) systematics due to problems with galaxy catalogues or d) 
intrinsic differences in the underlying galaxy clustering. In 
order to directly compare the 2dFGRS and SDSS, we anal- 
yse each dataset using essentially identical methods which 
we outline in Section [2] In Section [3] we briefly look at the 
region of overlap between the surveys to study the different 
selections used and the level of incompleteness. In Section 21 
we compare the resulting power spectra from the full cat- 
alogues and interpret the differences. Section [5] presents a 
detailed analysis of the implications of the differences be- 
tween these datasets in the obtained constraints on cosmo- 
logical parameters, and the systematic effects that might be 
introduced by the analysis. Finally, we conclude in Section[U] 



2 METHODS 

2.1 2dFGRS analysis 

The 2dFGRS covers approximately 1800 square degrees dis- 
tributed between two broad strips, one across the South 
Galactic Pole (SGP) and the other close to the North Galac- 
tic Pole (NGP), plus a set of 99 random 2 degree fields 
spread over the full southern galactic cap. The final cata- 
logue contains reliable redshifts for 221 414 galaxies selected 
to an extinct ion-corrected magn i tude l imit of approximately 
bj = 19.45 (|Colless et all 1200 ll I2003T I. The 2dFGRS sam- 
ples analysed he re are the same as described in detail in 
ICole et all (120051 ). 

Our method of estimating the galaxy power spectrum 
and determining s tatistical e r rors is essentially identical 
to that set out in ICole et al l l|2005l ). but with two minor 
changes. In brief: 

We use masks, whose construction is described in 
iNorberg et all (|2002T l. to describe the angular variation of 
the survey magnitude limit, redshift completeness and mag- 
nitude dependence of the redshift completeness. 

Random catalogues are generated by sampling from the 
luminosity function and viewing through the masks. To gen- 
erate random catalogues corresponding to red/blue subsets, 
the luminosity function of only the red/blue galaxies is used. 

Close pair incompleteness due to "fibre collisions" is 
dealt with by redistributing the weights of missed galaxies 
to their 10 nearest neighbours on the sky. 

The power spectrum is estimated using a simple cu- 
bic FFT method with the opti mal weighting scheme of 
iPercival Verde fc Peacock] l|2004 hereafter PVP) and then 
spheri cally averaged in redshift space. As in I Cole et all 
( |2005i ), the assumed linear empirical bias factors that are 
used in this weighting scheme are 

febiuc = 0.9 (0.85 + 0.15 L/L*), (1) 
for rest frame bj — r-p < 1.07, and 

6 rcd = 1.3 (0.85 + 0.15 L/L,), (2) 
for rest frame bj — tf > 1.07. These are based on the 



scale-independent bias parameter found by Norberg et al. 
(2001), modified to describe the difference in amplitudes 
of the power spectra of red and blue galaxies around k = 
0.1/iMpc - . This simple description of the luminosity and 
colour de pendence of ga l axy b ias is in agreement with the 
resu lts of ISwanson et all l|2007h . 

I Cole et al. I (|2005h showed that the effect on the recov- 
ered p ower spectrum of using iFeldman. Kaiser fc Peacock! 
(| 19941 hereafter FKP) rather than the PVP estimator is 
small (see their figures 17(s) and (t)). The FKP estimator 
is biased, as it ignores the luminosity dependence of galaxy 
clustering. Provided the model of luminosity-dependent bias 
is correct, then the PVP estimator removes this bias. A more 
complicated recipe for the luminosity and colour dependence 
of galaxy bias can be applied but, as long as the assumed bias 
is scale independent (a necessary assumption in our method 
of analysis) it will hardly change the shape of the estimated 
power spectrum, as the effect will be significantly smaller 
than the difference between the FKP and PVP estimates 
and so entirely negligible. 

The covariance matrix describing the errors on the 
power spectrum measurements and their correlations is es- 
timated using mock catalogues which are constructed from 
the random catalogues by generating a log-normal density 
field with a specified power spectrum and using it to mod- 
ulate the selection of galaxies from the random catalogue. 
Thus, by construction, these catalogues have a power spec- 
trum very close to the best fitting model and have luminos- 
ity and colour dependent clustering consistent with the bias 
factors of equations ([1} and ©. 

The survey window function is determined directly from 
the random catalogue. When fitting models the theoretical 
model power spectra are convolved with the survey window 
function. 

The two minor changes we have made are: 

• We have adopted a simpler binning scheme so that now 
P(k) is estimated in bins uniformly spaced in log 10 k, rather 
than the linearly space bins with dif ferent bin widths in dif- 
ferent ranges of k that were used in ICole et all l|2005l ). 

• We used new sets of log-normal catalogues in which the 
modulation of the density field used for galaxies with differ- 
ent bias factors is linear rather than t he slightly more co m- 
plicated scheme that was employed by I Cole et ail (|2005h . 



2.2 SDSS analysis 

We have analysed the SDSS Data Release 5 (DR5) sam- 
ple (|Adelman- McCarthy et aTll2007h . which is considerably 
larger than the DR3 sample analysed by iTegmark et all 
(2004). In most respects our analysis of the SDSS is identical 
to that of the 2dFGRS. The only differences are a simpler 
way of generating the survey masks and of populating the 
corresponding random catalogue. 

The sky coverage mask we have adopted for the SDSS- 
DR5 data is shown in Fig.[T]and compared with correspond- 
ing 2dFGRS mask. We constructed this mask by simply not- 
ing the angular coverage of each of the stripes from which 
the SDSS survey is built and by removing a few small re- 
gions with poor coverage. Most of the SDSS survey goes to 
a uniform magnitude limit of r = 17.77, but a sub-area, 
which is easily identified using the target selection date, has 
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Figure 2. The redshift distribution of the SDSS r < 17.77 galaxy 
sample (solid line). The dashed line shows an analytic fit to the 
data. The dot-dashed line is this same fit scaled down in ampli- 
tude to be always below the redshift histogram. This fit is used to 
generate the corresponding random catalogue by sampling from 
the original survey as described in the text. 



a variety of different magnitude limits. In this sub-area we 
simply imposed a fixed magnitude limit of r = 17.5 and 
discarded all galaxies fainter than this limit. The number 
of galaxies with redshifts that are retained by the mask 
and magnitude limits and used in our analysis is 443 424. 
The mask is cruder than t he mo re sophisticated ones em- 
ployed by iTegmark et~ai1 \2004 ) and for the 2dFGRS as 
it ignores the smaller scale variation in the redshift com- 
pleteness. However in the case of the 2 dFGRS, wh e re the 
incompleteness variation is much larger, ICole et al.l (I2005T I 
showed that provided this incompleteness is accounted for 
using our method of redistributing the weights of galaxies 
without redshifts to neighbours with redshifts, the result- 
ing power spectrum estimates are very accurate (see their 
figure 17g). 

As most of the catalogue goes uniformly to the deeper 
magnitude limit a simple method can be used to construct 
the random catalogue. Fig. [5] shows the redshift distribu- 
tion of this sample (solid line) together with an analytic fit 
(dashed line) that smoothes away the effect of large scale 
clustering. Our procedure is: 

(i) Select a random direction on the sky. 

(ii) Choose at random a genuine galaxy from the region 
of the survey that goes to r = 17.77 . 

(iii) Keep the galaxy with a probability proportional to 
the ratio of the fit scaled down in amplitude to be always be- 
low the redshift histogram (shown as the dot-dashed line in 
Fig. [2]| to that of the actual height of the redshift histogram 
at the redshift of the selected galaxy. This effectively adjusts 
the sampling of the different redshifts to erase the effect of 
large scale clustering. 



(iv) Keep or discard the galaxy based on the sky coverage 
and magnitude limit of the mask. (Note, for galaxies that fall 
where the magnitude limit is only 17.5, the fainter galaxies 
will be discarded and the redshift distribution of the retained 
galaxies will be appropriately shallower than that of Fig. [2]) 

These steps are done repeatedly until a random cata- 
logue containing 100 times more galaxies than the genuine 
catalogue is built up. Selecting from the genuine catalogue in 
this way means we automatically have apparent magnitudes 
and colours for all the galaxies in the random catalogue and 
so can select sub-samples from it and weight its galaxies in 
just the same way as the genuine catalogue. 

To utilize the iPercival. Verde fc Peacockl |2004l) opti- 
mum weighting we need to determine bias factors for the 
galaxies in the genuine, random and mock catalogues. We 
do this by first converting the SDSS magnitudes to the 2dF- 
GRS bj and tf bands using 



6j = g + 0.15 + 0.13 (g-r) 
r F = r-0.13 



(3) 



and the simple colour dependent fc-corr ections that wer e 
used for the 2dFGRS data (see section 3 of lCole et ai]|2005f ). 
Then we are able to define the bias factors using equa- 
tions JTJ and © just as for the 2dFGRS data. 



3 SURVEY OVERLAP 

Our main focus is the comparison of the power spectra of the 
two surveys, but we first directly compare the two surveys 
in the region of their overlap to get a feel for the different 
selections used and the level of incompleteness. 

In the northern galactic hemisphere there is a contigu- 
ous area of overlap between the two surveys, which runs for 
74 degrees of RA and is for the most part 5.2 degrees wide 
in declination. If we select from the SDSS photometric cat- 
alogue all galaxies brighter than bj = 20 (we do not apply 
the r « 17.77 magnitude limit of the SDSS main galaxy sur- 
vey, but we do apply all the oth er star-galaxy classi fication 
criteria used in that sample (see lStrauss et al. 2002)), then 
in this area there are 53 382 galaxies that are in both cata- 
logues. We find the fraction of SDSS galaxies which are also 
in the 2dFGRS to be constant at 89% as faint as fej R3 18.9. 
Fainter than this, SDSS galaxies are missing from the 2dF- 
GRS sample simply due to the (variable) magnitude limit 
of the 2dFGRS survey and its 0.15 magnitude random pho- 
tometric errors. This finding is in p erfect accord with the 
es timates made with th e SDSS EDR l|Stoughton et al.ll2002f ) 
in lNorberg et all (I2002T). Th i s 11% incomplet eness has been 
invest igated by Cross et all |2004h as well as iNorberg et al.l 
(2002) and has been shown to be predominately due to in- 
correct star-galaxy classification. The star-galaxy classifica- 
tion parameters based on the APM photometry are noisy 
and this level of incompleten ess is in line with what was 
expected l|Maddox et al.ll 19901 ). 

The issue here is whether this incompleteness has any 
influence on estimates of galaxy clustering. We can look at 
this directly by plotting cone plots (Fig. [3j of the galaxies 
the two catalogues have in common and those missed by 
the 2dFGRS. We plot only galaxies with bj < 19 to avoid 
issues with the 2dFGRS magnitude limit. We see that 91% 
of the SDSS galaxies are in the 2dFGRS. Comparing the two 
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Figure 3. Cone plots showing RA and redshift for galaxies in the region of sky where the 2dFGRS and SDSS surveys overlap. There 
are 53 382 within this area that are in both surveys and in the upper panel we plot the sub-sample of 27 004 that (estimated from SDSS 
photometry) have fej < 19 and redshift z > 0.01. The lower panel shows the 2754 galaxies that are in SDSS and pass the same magnitude 
and redshift cuts, but are missing from the 2dFGRS catalogue. 



cone plots in Fig. [31 it appears that the galaxies missed by 
2dFGRS are just a random sparse sampling of the structure 
seen in the matching sample and so there is no evidence that 
the incompleteness has a spatial imprint. 



4 COMPARISON OF POWER SPECTRA 

The power spectra we estimate using the methods outlined 
in Section[5]are the direct transform of the data, and are thus 
what CMB researchers would term a pseudo-spectrum. As 
such, it yields a convolution of the underlying galaxy power 



spectrum with the modulus squared of the Fourier transform 
of the window function of either the SDSS or 2dFGRS as 
appropriate, that is 

P(k) = P(k) <g> W 2 (k). (4) 

Our random catalogues allow us to accurately estimate 
VF 2 (k) and from it determine the matrix of window func- 
tions that describe how our spherically averaged band power 
estimates P(kj) are related to the unconvolved power spec- 
trum P(ki). 

P(kj) = ^2P(h)W m (h,k j ). (5) 
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Figure 4. A sample of the normalized window functions for 
both SDSS (dashed) and 2dFGRS (solid). Each curve shows the 
relative contribution from the underlying power spectrum P(k{) 
as a function of ki to our individual band power estimates P(kj). 
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Figure 5. Comparison of the power spectra estimated from the 
full 2dFGRS and SDSS DR5 samples - corrected for the effect 
of the window function as described in the text. The solid line 
shows the input power spectrum of the mock catalogues used to 
estimate the covariance matrix of the measurements. 



Examples of these window functions for the SDSS and 2dF- 
GRS are shown in Fig. [4] For all our quantitative analysis 
in Section [S]we use these window functions to convolve the 
model power spectra before comparing with the data. How- 
ever for the purposes of visually comparing the 2dFGRS and 
SDSS power spectra we have corrected the convolved esti- 
mates by multiplying them through by the ratio of a similar 
model power spectrum to its convolved counterpart. This 
'deconvolution' is accurate provided the power spectra are 
smooth. 

In Fig. [S]we compare the 'deconvolved' power spectra 
estimated from the full 2dFGRS catalogue and full SDSS- 
DR5 sample. We note that the SDSS and 2dFGRS galaxy 
power spectra agree well for k > 0.07 ft Mpc -1 . The good 
agreement in amplitude at this wavenumber is due to the 
bias dependent weights used in the PVP estimator which 
have successfully modelled the difference in the clustering 
strength of the red selected SDSS galaxies and blue selected 
2dFGRS galaxies. This is by design as the bias factors were 
normalized empirically by the 2dFGRS red and blue samples 
at this scale (see figure 15 of ICole et al]|2005h . In contrast, 
on larger scales we see evidence for significantly more large 
scale power in the 2dFGRS than in SDSS. This clearly shows 
that the difference between the SDSS and 2dFGRS results 
that was noted in the introduction and which motivated 
this analysis is certainly significant and not an artifact of 
differing analysis techniques. 

We now investigate if the discrepant shapes of the 
galaxy power spectra are due to the difference in the clus- 
tering properties of red and blue galaxies. Fig. [6] shows his- 
tograms of rest frame bj — tf colours for both the 2dFGRS 
and SDSS catalogues. The SDSS magnitudes have been con- 
verted to these bands assuming the relations given in equa- 
tion ((3|. The colour distributions are clearly bimodal with 



a natural dividing point at bj — tf = 1.07. The 2dFGRS 
has roughly equal numbers of red and blue galaxies while 
the SDSS, being red selected, is naturally dominated by red 
galaxies. 

In Fig. [7] we compare the galaxy power spectra esti- 
mated from just the galaxies redder than bj — rp = 1.07 
in both samples. Comparing the SDSS P(k) from just the 
red galaxies with the previous estimate from the full SDSS 
catalogue reveals them to be in very close agreement. This 
is to be expected as the SDSS sample is both dominated by 
red galaxies and the PVP power spectrum estimator that 
we employ gives them more weight than their less clustered 
blue counterparts. In contrast, the estimate for just the red 
2dFGRS galaxies differs from that from all the 2dFGRS cat- 
alogue. In fact it is a much closer match to the result from 
the SDSS data. The only places where the two estimates 
are not in excellent agreement is on the very largest scales 
k < 0.025 h Mpc -1 , where the estimates are both noisy and 
highly correlated, and also around k w 0.05 h Mpc -1 . In 
fact th is difference is also due to sample variance. I Cole et all 
(2005) investigated the effect of removing from the 2dFGRS 
catalo gue the two l argest super clusters found by the anal- 
ysis of iBaugh et al.l l|2004l ). Their figure 17 (panels o and p) 
shows that this in general has a small effect, but does per- 
turb the power just around k « 0.05 /iMpc" . 



5 IMPLICATIONS FOR COSMOLOGICAL 
PARAMETERS 

5.1 The shape of the power spectrum 

The power spectrum measured for galaxies differs in a num- 
ber of ways from the power spectrum for the mass predicted 
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Figure 6. Histograms comparing the distribution of rest frame 
6j— rp colours in the 2dFGRS (solid line) and SDSS (dashed line) 
catalogues. The distributions are clearly bimodal with a dividing 
point at fej — rp = 1.07. The 2dFGRS has roughly equal numbers 
of red and blue galaxies while the SDSS is dominated by red 
galaxies. Note the different units on the y-axis. 



Figure 7. Comparison of the 2dFGRS and SDSS-DR5 power 
spectra from red subsamples that satisfy the rest frame colour 
6j — rp > 1.07. The solid line shows the input power spectrum 
of the mock catalogues used to estimate the covariance matrix of 
the measurements. 



in linear perturbation theory: (i) Nonlinear evolution of den- 
sity perturbations leads to coupling between Fourier modes, 
changing the shape of the power spectrum, (ii) The galaxy 
power spectrum is distorted by the gravitationally induced 
peculiar motions of galaxies when a redshift is used to infer 
the distance to each galaxy, (iii) The power spectrum of the 
galaxies is a modified version of the power spectrum of the 
mass. This phenomenon is known as galaxy bias. As we have 
shown in the previous section these effects change with scale 
and depend on the galaxy type. 

In order to constrain c osmological param eters, these ef- 
fects need to be modelled. I Cole et all (|2005t ) developed an 
empirical scheme to correct for these effects by applying a 
correction for non-linearity and scale-dependent bias to the 
shape of P(k) of the form 

P g Ak) = b 2 ^±^P lin (k), (6) 

where Pu n (k) is the linear theory power spectrum and b is a 
constant bias factor. This formula was deduced by compari- 
son with detailed numerical galaxy-formation models: these 
show that the value of A = 1.4 is robust, but the exact value 
of Q depends on galaxy type and also has some uncertainty 
depending on how the modelling is done. These results were 
used to determine a range of allowed values for Q. For 2dF- 
GRS the value Q = 4.6 is preferred. 

Although this correction factor was designed and tested 
for the redshift space power spectrum of 2dFGRS where 
the scale dependent correction is small, it has also been 
used to model the real space p ower spectra for redder and 
more luminous galax y samples ijTeemark et al.ll2004l . [200(3 : 
IPercival et al. 2007a b). In these cases larger values of Q 
were required to reconcile the measured power spectrum 



with linear theory (e.g Q = 26 in IPercival et all l2007tj ). 
Here we analyse the validity of this correction in general 
terms, testing if the different shapes of the power spectra 
for the full and red galaxy samples of 2dFGRS and SDSS 
can be accounted for by using equation (0 with different 
values of Q. 



Besides the change in the overall shape of P(k), it has 
been shown that non-linear evolution distorts the acous- 



tic oscillations ( 


Eisenstein et al. 2005^ SDringel et al.ll2005l: 


lEisenstein et al. 


2006!; lAnsulo et al. 2005, 20071). These os- 



dilations are damped in a way that erases the higher har- 
monic peaks. lEisenstein et al.l (|2006l ) modelled the damping 
of the oscillations by computing a 'dewiggled' power spec- 
trum by 



Pdw(fc) = P lin (k)G(k) + P nw (fc)(l - G(k)), (7) 

where P nv ,(k) is a smooth power spectrum, with the same 
shape as Py ln (k) but without b aryonic oscillation s , comp uted 
using the fitting formulas of lEisenstein fc Hvj (|l999l ) and 
G(k) = exp [— (fc/\/2fc*) 2 ] • This function regulates the tran- 
sition from the large scales, where Pdw(fc) follows linear the- 
ory to the small scales where the a coustic oscillations ar e 
completely damped. As described in lTegmark et ail {2006), 
the value of the damping scale fc* is a function of O m and 
the primordial amplitude of scalar fluctuations A s . Here we 
analyse the effect of this modelling on the obtained con- 
straints on cosmological parameters. When this correction 
is applied, Pdw(k) is used in equation instead of Pn n (k). 

We use the data from the observed power spectra for 
k < 0.2/iMpc" 1 and discard measurements with k < 
0.02 ZiMpc" 1 which could be affected by uncertainties in the 
mean density of galaxies within the survey. We compare the 
data against a restricted parameter space given by 
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Figure 8. Contours show 68% and 95% joint confidence intervals 
for the baryon fraction f2 b /f2 m and f2 m /i for fits to the SDSS 
and 2dFGRS data in the range 0.02 < k < 0.2/iMpc -1 . The 
parameter Q modelling the distortion of power spectrum due to 
nonlinearity, redshift space distortions and scale dependent bias 
was kept fixed at Q = 5 in these fits. 



p = (n m ft,n b /n m ,Q, J 4 8 ), (8) 

and we fix the values of h = 0.72 and n s = 1. Eq. ((5} gives 
the most important parameters to characterize the full shape 
of P(k). When constraining the values of these parameters, 
we only use information from the shape of the galaxy power 
spectrum and not from its amplitude. This is why A 3 only 
enters in the constraints marginally when the damping of the 
acoustic oscillations is included in the mode lling through the 
assumed dependency of fc* on fi m and A s (|Eisenstein et al.1 
l200fj ; lTegmark et al"j|20o3 ). 



Our results were generated with a modified version o f 
the publicly availab l e Cos moMC code (|Lewis fc Bridlell2002r ) 
as in lSanchez et all (120061 ). CosmoMC uses the CAMB pack- 
age to co mpute the linear p ower spectra for the matter fluc- 
tuations dLewis et al.ll2000[). We have fou nd that the use of 
the approximate lEisenstein fc Hul l|l998l ) fitting formula to 
model Pnn(fc) causes a small shift in the recovered values of 
f2 m /i. Our analysis was carried out in parallel on the Cos- 
mology Machine at Durham University. For each parameter 
set considered, we ran 10 separate chains using the Message 
Passing Interfac e (MPI) convergence cr iterion to stop the 
chains when the lGelman fc Rubin! l|l992h statistic R < 1.01, 
which is a sign ificantly more stringent criterion th an is usu- 
ally adopted (|Verde et al.ll2003l ; ISeliak et alJl200Sj ). 

To test the efficiency of the modelling to account for 
the differences between the datasets we analyse three sepa- 
rate cases: first we fixed the parameters A and Q at fiducial 
values of A = 1.4 and Q — 5 and ignore the correction due 
to the damping of the acoustic oscillations (Section 15. 2[) . In 
this way the difference in the shape of the power spectra for 
SDSS and 2dFGRS is completely translated into the recov- 
ered constraints on fi m /i and fib/fim. Second we extend the 
model allowing Q to vary over a wide prior (Section 15.31) . 
In this way we can test if the ansatz of equation © can 
account for the observed differences. Finally, we include the 
effect of the damping of the acoustic oscillations into our 
modelling of P(k) (Section I5.4[) . Table Q] summarizes the 
constraints obtained with the different datasets and param- 
eter sets analysed. 

5.2 Fixed Q 

First we analyse the constraints on Q m h and f2b/fi m that 
are obtained by fitting the SDSS and 2dFGRS data with 
power spectra of the form described in the last section, with 
the parameters A and Q fixed at fiducial values of A — 1.4 
and Q = 5 and ignoring the correction due to the damp- 
ing of the acoustic oscillations. In this way these are the 
only two parameters defining the shape of the model P(k) 
and the difference in shape of the fitted power spectra is 
completely characterized by the regions of the explored pa- 
rameter space selected by each dataset. The parameter con- 
straints obtained in this way are shown in Fig. [5] for the full 
(upper panel) and red (lower panel) 2dFGRS (solid lines) 
and SDSS (dashed lines) galaxy power spectra. 

Fig. Ufa) shows that the 2dFGRS constraints are in 
complete agreement with those obtain ed by CMB data 
|Sanchez et al.ll200d ; lsp"ergel et al.ll2007f l. For 2dFGRS data 
only w e obtain Q. m = 0.246 ±0.030 while the WMAP3 anal- 
ysis of ISpergel etal] (|2007f l gives Q, m = 0.234 ± 0.035. The 
SDSS para meter estimates are c ompletely in accord with 
those from iTegmark et al.l ([2004) , but with much tighter 
bounds due to the larger SDSS dataset used here. This shows 
the consistency of the different analysis techniques applied 
to this dataset. 

We note that the 2dFGRS and SDSS best fit values lie 
outside each others 95% confidence contours. This clearly 
shows that, if not accounted for, the difference between the 
SDSS and 2dFGRS power spectra that was noted in Sec- 
tion 0] causes a bias in the obtained constraints on cosmo- 
logical parameters which is larger than the uncertainty with 
which they are determined. 



Precision Cosmology from Large Scale Structure? 9 

Table 1. Marginalized 68% interval constraints on cosmological parameters obtained using the different datasets and parameter sets 
analysed. 



Data set 


Parameter 


Fixed Q 


Varying Q 


Varying Q, dw 


2dFGRS all 


Q 


n 177 +0.022 
u - 11 ' -0.025 

U ' i95 -0.057 


n lfiK+0.025 
u.ido_ 02g 

n 91 9 +0.061 
u - ziz -0.065 
7 7+3.3 
'■'-3.4 


191+0-064 
U.±M1_ q 51 

1+3.7 
8 " L -3.8 


2dFGRS red 


Q 


n 91 7 +0.040 
U - Zi ' -0.040 
n 1QQ+0 064 


iqy+0.043 
u - ±31 -0.043 

0.197™ 


023+ 094 

fl 9Q0+O.O86 
u.zoo_ Q Qg9 


SDSS all 


Cl m h 

Q 


199+ ' 014 

n 19H +0.029 
U-128_ 031 


216+ 019 

U.Z1D_ 02Q 

o.nai° :S§ 


0.234±° : ™ 

n i= n +0.050 

u.iou_ 053 

1+1.6 


SDSS red 


f^b/^m 

Q 


201+ 013 
U ' ZUi -0.015 

120+ 027 
u.i^u_ 030 


221+ 017 
U - ZZi -0.025 

i ,-+0.031 

U - liO -0.036 

o 9+I.6 
cs.z_ 1 7 


240+ 038 
u.^u_ Q 039 

n -,,-7+0.057 
U - iO '-0.058 

3.5 17 


LRG 


^b /£lm 

Q 




n i 79+O.OI6 
u ' 1( -0.016 
n i 7 o+0. 038 
ul '°-0.036 

27 1+ 4 ' 7 


Q -,73+0.017 

u - x '°-0.017 

0.187+°;^ 
28.51" 



Fig. [8fb) confirms the agreement, within the expected 
statistical uncertainty, of the power spectra for red galaxies 
in 2dFGRS and SDSS. Being dominated by red galaxies, the 
SDSS constraints are almost identical to the ones obtained 
with the full sample. Instead, the results for 2dFGRS shift 
towards the same region of the parameter space preferred by 
SDSS to contain it completely. The best fitting parameters 
for 2dFGRS lie within the SDSS 65% confidence contour and 
vice versa. The contours obtained for 2dFGRS broaden as 
a result of the larger errors in the power spectrum obtained 
with a smaller sample of galaxies. This shows that the differ- 
ences in the constraints obtained with SDSS and 2dFGRS 
are due to the different selection criteria used to construct 
the samples, with the r-band selected SDSS catalogue being 
dominated by more strongly clustered red galaxies which are 
more strongly affected by scale dependent bias. 



5.3 Varying Q 

In this section we extend the parameter space by allowing 
Q to vary over a wide prior. This allows us to test if the 
shape distortion modelled by equation (JSJ can account for 
the differences between 2dFGRS and SDSS found in Sec- 
tion U The results obtained in this way are shown in Fig. 
for the full (upper panels) and red (lower panels) 2dFGRS 
(solid lines) and SDSS (dashed lines) galaxy power spectra. 
For this case, we have also included in our analysis the power 
spectrum of the SDSS Luminou s Red Galaxies (LRG, dotted 
lines in Fig. [9jl as measured bv lTegmark et all |2006). 

From Figs. |9ja) and (b) we note that the correction fac- 



tor of equation © is not able to reconcile the constraints 
obtained from the 2dFGRS and SDSS power spectra. The 
constraints on the STi m — fib/fim plane exhibit a similar be- 
haviour to the case where the value of Q was held fixed. This 
shows that the non-linearities and scale dependent bias of 
the r-selected SDSS galaxies distort the shape of the power 
spectrum in a way that can not be described accurately by 
equation ©. Contrary to what was expected, the constraints 
on the plane f2 m /i — Q show that SDSS and 2dFGRS select 
different values of Qmh but a similar allowed region for Q. 
Then, even incorporating this model of the distortion into 
the methodology for constraining cosmological parameters, 
the use of a red sample of galaxies introduces strong system- 
atics effects on the resulting parameter constraints. 

Intriguingly, the constraints on Q m h and fib/f^m ob- 
tained using the LRG power spectrum lie well within the 
2dFGRS contours, in perfect agreement with the results ob- 
tained using CMB data. This means that the difference in 
the shape of the SDSS-LRG and 2dFGRS power spectra 
is completely consistent with equation (HJJ. Both datasets 
can be accurately described with the same values of the 
cosmological parameters, but with a much higher value of 
Q = 27.1±4. 7 for the SDSS-LRG sa mple, in agreement with 
the results of iPercival et all (|2007bh . 

As shown in Fig. H^c) and (d), the situation is similar 
for the constraints obtained using the power spectra of the 
red 2dFGRS and SDSS galaxy samples. As in the case of a 
fixed value of Q, the constraints obtained with the 2dFGRS 
power spectrum are shifted from its region of agreement with 
the CMB towards the values preferred by the SDSS. With 
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Figure 9. Contours show 68% and 95% joint confidence intervals for the baryon fraction £7]->/f2 m and £7 m /i for fits to the 2dFGRS (solid 
line), SDSS (dashed line) and SDSS-LRG (dotted line) in the range 0.02 < k < 0.2 h Mpc~ x . In these fits, the parameter Q modelling 
the distortion of power spectrum due to nonlinearity, redshift space distortions and scale dependent bias was allowed to vary over a wide 
prior (0 - 45). 



the addition of Q, the extra degree of freedom gives rise to 
a degeneracy that involves fi m /i and Q. These two parame- 
ters act in opposite ways which allows them to compensate 
each others effect on the theoretical power spectrum. This 
degeneracy can not be broken by the power spectrum of the 
red 2dFGRS galaxies and results in a considerable increase 
of the allowed region for the different parameters. 

The breakdown of the simple model of equation (JSJ) in 
describing the differences in the clustering of red and blue 
galaxies can be more clearly seen in Fig. 1101 which shows 



the different power spectra analysed in this work divided 
by a reference power spectrum P Te {(k ) computed us i ng th e 
best fit cosmological parameters from lSpergel et all l|2007l) . 
Assuming this is the true cosmology, if there were no dis- 
tortions in the shape of P(k), the different measurements 
should lie over the dashed line of P(k)/P IC ; (k) = 1. The dif- 
ference in the shapes of these spectra is a clear demonstra- 
tion of the fact that scale dependent bias and nonlinearities 
are a function of galaxy type. The solid lines in Fig. HOI show 
the shape distortion model of equation ((6]) for the specified 
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Figure 10. Th e power spectra anal ysed in this work divided by a reference power spectrum computed using the best fit cosmological 
parameters from lSpergel et alj d2007t) ■ The solid lines show the shape distortion modelled by equation (JBJl for the specified value of Q. 



value of Q. With a value of Q ~ 5 this correction gives a 
good description of the distortions in the power spectrum 
of the 2dFGRS, with a small decrease in the amplitude of 
P(k) at intermediate scales and an increase at small scales. 
The same correction factor can reproduce the effect of the 
larger non-linearities of the LRG power spectrum, but with 
a much higher value of Q ~ 30. 

Despite its success at describing the 2dFGRS and LRG 



power spectra, equation fails to reproduce the shape of 
the SDSS power spectrum, which exhibits a different behav- 
ior as a function of k and can not be fitted by this simple 
model. The fact that the power spectrum of the red galaxies 
from 2dFGRS is distorted in a similar fashion indicates that, 
as was found in Section |4j red selected galaxy samples ex- 
hibit stronger scale dependent galaxy bias. For comparison, 
Fig. [10] also shows that the power spectrum of SDSS DR3 
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measured by iTegmark et"al] (|2004T l suffers from the same 
effect. 

This clearly shows that, in order to obtain unbiased con- 
straints on cosmological parameters it is necessary to better 
understand and constrain the distortion in the shape of the 
power spectrum caused by nonlinearity and scale dependent 
bias. 



5.4 Damping of acoustic oscillations 

In this section we analyse the effect of including the damp- 
ing of the acoustic oscillations in to our modelling of P (k). 
This case mimics the analysis of ITegmark et all (|2006T l of 
the LRG power spectrum. Although we only use informa- 
tion from the shape of the galaxy power spectrum, when 
this effect is included in our modelling of P(k) the am- 
plitude As enters in the constraints marginally, since the 
value of the damping scale fc* depends on this parameter 
l|Eisenstein et al.ll2006l ; ITegmark et alj|2006h . 

As can be seen in the last column of Table Q] this addi- 
tion does not change the differences between 2dFGRS and 
SDSS. The main effect of this is the increase of the allowed 
regions of the different cosmological parameters. As the os- 
cillations are damped, models with higher values of fib/f2m, 
which according to linear theory exhibit larger amplitude 
oscillations, are now compatible with the data. This causes 
the signal of the oscillations to effectively loose the power 
to constrain efficiently the values of the parameters, which 
in this case are determined almost completely by the overall 
shape of P(k). 

The details of the way in which the oscillations 
are damped and the correct way to model it are the 
main focus of several recent studie s (Eisen stein et al]|2006l ; 
lAngulo et aHl2007l ; ICrocce fc Scoccimarrol 120071 ). Although 
substantial progress has been made in this subject, it is not 
completely understood. This is then another example of how 
the details of the modelling used to fit theoretical models 
may introduce systematic effects which may cause signifi- 
cant differences in the obtained constraints. 



6 CONCLUSIONS 

In this paper we studied in detail the differences between 
the constraints on cosmological parameters coming from 
the 2dFGRS and SDSS power spectra. We have analysed 
both the 2dFGRS and SDSS DR5 catalogues using essen- 
tially identical techniques. The main conclusion of this in- 
vestigation is that a significant difference exists between the 
shape of the galaxy power spectra measured in the 2dF- 
GRS and SDSS main galaxy sample and that this difference 
is due to scale depende nt bias. Our results extend those 
of IPercival et al.l |2007a) who found that the inclusion of 
the LRGs in the analysis increases the discrepancy between 
2dFGRS and SDSS. Their findings also showed that if lin- 
ear CDM fits are restricted to large scales (0.01 < k < 
0.06 h Mpc^ 1 ), SDSS dat a also favours a lower m atter den- 
sity Q, m h = 0.16 ± 0.03. IPercival et all l|2007al ) suggested 
that a scale-dependent bias which depends on the r-band 
luminosity of the galaxies could explain these differences, 
but found no significant evidence of it. 

If a homogeneous sample of red galaxies is selected from 



each survey then the resulting power spectra agree to within 
the expected statistical errors. In contrast, when the full 
2dFGRS and SDSS catalogues are analyzed the resulting 
2dFGRS power spectrum differs in shape to that of SDSS. 
If normalized on scales around k w 0.1 ftMpc" 1 , as is done 
automatically by the scale independent bias factors assumed 
in the PVP estimator, then the 2dFGRS P(k) exhibits more 
large scale power than SDSS. However, if one instead nor- 
malizes on large scales one finds the equivalent result that 
SDSS exhibits more small scale power than 2dFGRS. This 
behaviour is exactly what one expects if the more strongly 
clustered red galaxies live in denser environments where the 
effects of nonlinearity are greater. 

We have obtained constraints on the parameter com- 
binations Qmh and Qb/Q m which are the most important 
parameters to characterize the full shape of P(k). The con- 
straints obtained by fitting the SDSS and 2dFGRS data 
show that the difference between the SDSS and 2dFGRS re- 
sults which motivated this analysis is significant and is not 
due to the differing analysis techniques. If not accounted for, 
the difference between these datasets introduces a system- 
atic effect that bias the obtained constraints on cosmological 
parameters, whose effect is larger than the uncertainty with 
whi ch they are determ ined. 

I Cole et al] (|2005h introduced a first attempt at mod- 
elling the distortion in the shape of the power spectrum 
caused by nonlinearity and scale dependent bias using the 
Q and A parameters of equation For the mix of red 
and blue galaxies present in the 2dFGRS survey, the nec- 
essary value of Q was reasonably small and hence the scale 
dependent correction quite modest. The same correction fac- 
tor has been applied to samples of redder or more luminous 
and hence more clustered galaxies where one expects greater 
nonlinearity. 

We found that the correction factor of equation ([6]) is 
not able to reconcile the constraints obtained from the 2dF- 
GRS and SDSS power spectra. This shows that this simple 
correction can not describe the effect of non-linearities and 
scale dependent bias for a general sample of galaxies. In- 
triguingly, the same model can describe correctly the differ- 
ence in the shape of the power spectra from 2dFGRS and the 
much strongly clustered SDSS-LRG sample. Both datasets 
can be described with the same values of the cosmological 
parameters, but with a much higher value of Q = 27.1 ± 4.7 
for the SDSS-LRG sample. 

This comparison has revealed that to get unbiased es- 
timates of the cosmological parameters it is necessary to 
better understand and constrain the different processes that 
shape the galaxy power spectrum. Thus to get robust con- 
straints from the main SDSS survey and from future galaxy 
surveys like Pan-STARRS, or the Dark Energy Survey, 
which will use selection criteria similar to the one of SDSS, 
will require more detailed modelling of nonlinearity and scale 
dependent bias. 
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